Signature of quantum criticality in photo-emission spectroscopy 
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A quantum phase transition (QPT) in a heavy-fermion (HF) compound may destroy the Fermi 
liquid groundstate. However, the conditions for this breakdown have remained obscure. We report 
the first direct investigation of heavy quasiparticle formation and breakdown in the canonical sys- 
tem CeCue-^Auj; by ultraviolet photoemission spectroscopy at elevated temperatures without the 
complications of lattice coherence. Surprisingly, the single-ion Kondo energy scale Tk exhibits an 
abrupt step near the quantum critical Au concentration of x c = 0.1. We show theoretically that this 
step is expected from a highly non-linear renormalization of the local spin coupling at each Ce site, 
induced by spin fluctuations on neighboring sites. It provides a general high-temperature indicator 
for HF quasiparticle breakdown at a QPT. 
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PACS numbers: 71.27.+a 71.28. +d 79.60.-i 71.10.-w 

A quantum phase transition (QPT) is a second-order 
transition occurring at zero temperature, driven by a 
nonthermal control parameter such as composition or 
pressure. Here, two competing groundstates are sepa- 
rated by a quantum critical point (QCP). Heavy-fermion 
(HF) systems, in particular, show striking deviations 
from Fermi-liquid (FL) behavior at a QPC between a 
magnetically ordered and a paramagnetic phase 

[Ell- In 

the standard Hertz-Millis (HM) scenario 0, H only the 
bosonic fluctuations of the order parameter become quan- 
tum critical, i.e., long-ranged in space and time, leading 
to anomalous behavior of physical quantities, but leaving 
the fermionic quasiparticles intact. However, in a heavy 
fermion (HF) system the order-parameter fluctuations al- 
ways couple to the spin degree of freedom of the fermionic 
excitations, so that the latter may become critical as well 
and disintegrate. In this case, the formation of the Kondo 
spin singlet between the conduction electron and the lo- 
cal 4/ magnetic moments is prevented, and hence the 
Kondo scale Tk vanishes at the quantum critical point 
(QCP). This scenario has, therefore, been termed local 
quantum critical (LQC) @,@| CeCuR ^Au x is one of the 
best characterized HF compounds @, 0, El EH, El El 
with a QPT between a fully Kondo-screened paramag- 
netic and an antiferromagnetically ordered phase at a 
critical Au concentration of x c = 0.1. Based on inelastic 
neutron scattering experiments it was suggested that the 
LQC scenario is realized in this compound [3], yet this 
issue has remained controversial. 

In order to address the crucial question whether Kondo 
screening persists or not at the QPT, we present in this 
work ultraviolet photoemission spectroscopy (UPS) mea- 
surements, which provide the most immediate access to 
the screening scale Tk by directly recording the Kondo 
resonance (KR) in the local Cc 4/ spectrum [HI, flij |. 



By performing the experiments at temperatures T well 
above the Need temperature and above the lattice coher- 
ence temperature (T > {Tk, T co h] Tjy}) we are able to 
investigate the local spin screening effect alone, without 
the complications of Fermi-volume change due to lattice 
coherence, which influence the results of other experi- 
mental techniques. 

The high-resolution UPS was performed with a Gam- 
madata R4000 analyzer and a monochromatized VUV- 
lamp at hv = 40.8 eV (see Refs. [II and [13] for more 
details). We cleaved the single crystalline samples in situ 
just before the measurement, already at the measure- 
ment temperature. The surface is known to be y-Ce-like, 
meaning that the surface properties do not differ from 
the bulk properties [l6T |, 

The inset of Fig. QJi) shows survey spectra for five Au 
concentrations x at T=15 K. The energy range includes 
the spin-orbit (SO) feature at Eb ~ 260 meV and the 
tail of the Kondo resonance just below the Fermi energy 
Ep. A distinct crystal-field (CF) feature as observed, 
e.g., in CeCu2Si2 and CeNi 2 Ge2 [lallal. is expected near 
Eb = 7 meV, but strongly broadened [12| and, hence, not 
discernible in CeCug-a; Au^. The SO peak does not show 
a significant x-dependent variation in position or inten- 
sity. This is expected, since it corresponds to an only vir- 
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tually excited initial state of the UPS process 15, 18]. In 
order to investigate the KR in more detail we performed 
high- resolution measurements near Ep (Fig. Wfa)). For 
all concentrations one observes a significant rise in the 
intensity from Eb ~ 20 meV towards Ep, attributed to 
the tail of the KR which has its major spectral weight 
a few meV above the Fermi energy [la . 1 1 61 ] suppressed 
by the Fermi-Dirac distribution function (FDD) in the 
UPS data. Normalizing to the FDD allows to recover 
the thermally occupied spectrum up to sa 5ksT above 
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FIG. 1: a) Near-_E.F spectra of CeCu6_ x Aua; for five different Au concentrations at T—15 K (hv = 40.8 eVn AE — 4.9 meV), 
normalized at Eb ~ 100 meV. The dashed lines describe the resolution broadened FDD at T— 15 K. The inset shows a larger 
energy range including the spin-orbit (SO) partner at _Eb~260 meV (AU=15 meV)). b) and c) show spectra for £ = 0.1 and 
x = 0.2, respectively, devided by the FDD, at various T. The solid lines are best NCA fits (see Tab. [T] for the resulting model 
parameter values). The insets in b) and c) show the corresponding raw data. 



E F [3 Ell and reveals the KR. Fig. Q] shows the FDD- 
renormalized spectra for x = 0.1 (b)) and 0.2 (c)). In 
the low-T spectra (T = 15 K) its maximum is at about 
3 meV above Ep for x = 0.1 and at 1 meV for x = 0.2. 
The striking feature seen in Figs.[T]b) and c) is the signif- 
icant drop of the spectral weight of the KR from x = 0.1 
to 0.2, i.e., in the vicinity of the critical concentration x c . 
The differences in the shape and the T dependence be- 
tween the x = and 0.1 spectra as well as the differences 
between the x = 0.2, 0.5, and 1.0 spectra (not shown) are 
negligibly small compared to the jump close to x c . 

To extract the Kondo screening scale Tk from the ex- 
perimental results, one must bear in mind that the coher- 
ence temperature T co ^ L obtained from resistivity measure- 
ments [19B as well as the Kondo scale extracted from spe- 
cific heat (20l | or neutron scattering data II, U | , are well 
below our lowest experimental temperature. Moreover, 
quantum critical fluctuations, which in CeCu6-£cAu x ex- 
tend up to about T = 7 K [lj] and certainly become 
crucial at low T, are not expected to influence our exper- 
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TABLE I: SIAM/NCA fit parameters, CF splittings of the 
Ce 4/ orbitals, Acfi, Acf2, 4/-conduction band hybridiza- 
tion matrix element, V, and the resulting Kondo tempera- 
tures, Tk, for different Au concentrations x. Fixed param- 
eters: conduction-band half- width D = 2.8 eV (HWHM), 
single-particle 4/ binding energy e/ =— 1.05 eV, SO splitting 
A SO = 250 meV. 



imental spectra. Thus, our data exhibit the onset of the 
local Kondo physics on the Ce atoms only. Therefore, 
the single-impurity Anderson model (SIAM) is employed 
to interpret the experimental data. To determine Tk we 
follow the procedure successfully app lied to various Ce 
compounds in the past [H [H [2l| 



2^, HI: Using the 



noncrossing approximation (NCA) [25j we calculate the 
Ce 4/ spectral function of the SIAM, including all CF 
and SO excitations. For each composition x the NCA 
spectra are broadened by the experimental resulotion and 
fitted to the experimental data, using a single parameter 
set for all experimental T. The NCA spectra are then 
computationally extrapolated to TssO.l Tk, where Tk is 
extracted from the Kondo-peak half-width at half max- 
imum (HWHM) of the NCA spectra. The resulting fit 
parameter values are shown in Tab. [U Fig. [2] shows the 
x dependence of Tk for CeCug-zAu,; obtained from our 
UPS data, compared to results of various other experi- 
|. We emphasize that, irrespective of 
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a possible systematic ambiguity in the fit procedure, the 
surprisingly abrupt step of Tk near x = x c ~0.1 is signif- 
icant and already clearly visible in the raw data (Fig. [I}. 

What can be learnt from this step-like behavior at 
elevated T about the nature of the QPT at T = K? 
The method of extracting a Kondo scale from a single- 
ion picture described above mimics exactly a system of 
dense Kondo atoms that, however, would not undergo 
any lattice coherence or magnetic ordering at sufficiently 
low T. This scenario is different from the low-T behav- 
ior of the Kondo lattice and two-impurity Kondo prob- 
lems, where Kondo spin screnning would have to com- 
pete with the formation of inter-impurity coherent states 
2rj | or could be destroyed by critical fluctuations of a 
magnetic order parameter [27|. By contrast, our sys- 
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FIG. 2: Dependence of the Kondo temperature Tk on the Au 
content x, as determined by UPS (open circles), specific heat 
[20I ] (triangles) and neutron scattering [l2l H3 ] (diamonds). 
The shaded area is a guide to the eye. The insets a) and b) 
show diagrammatic representations of the RKKY corrections 
to the local Kondo vertex. The inset c) and the solid line 
in the main panel show the universal curve Tk(u)/Tk(0) vs. 
y/y m as given by Eq. Q. 

tem is represented by an effective single-impurity Kondo 
model where the Kondo exchange coupling J is renor- 
malized by the local spin fluctuations on the surrounding 
identical Kondo atoms through the indirect Ruderman- 
Kittel-Kasuya-Yosida (RKKY) coupling. This problem 
can be treated in a controlled way by renormalization 
group (RG) enhanced perturbation theory which allows 
for definite predictions about the formation of HF quasi- 
particles at low T, as we will show in the following. 

The leading perturbative corrections to J on a given Ce 
atom at lattice site i — are depicted diagrammatically 
in the insets a), b) of Fig. O They obviously involve the 
full dynamical spin susceptibility X4f(T, SI) of a local Ce 
4f spin on a neighboring site i ^ 0, where %4/(T, 0) — 

(9lLIb) 2 iV(0)ZV(V T A' + T2 )i with N (°) the density 
of states at the Fermi level, Dq « Ep the band cutoff 
and gL, the Lande factor and the Bohr magneton, 
respectively. Summing over all lattice sites i ^ 0, one 
obtains for the direct and the exchange corrections (Fig. [2] 
a) and b)), respectively, the spin-isotropic terms, 




SJ (ex) 



Here J is the bare, local spin exchange coupling on site 
i = 0, and gi = N(0)Ji is the dimensionless, bare coupling 



on site i ^ 0. y is a dimensionless factor that describes the 
(experimentally not known) relation between the RKKY 
coupling strength and the Au content x. We assume a 
linear dependence across the QPT, y = a(x + xq), with 
adjustable parameters a and xq. As seen from Fig. [3 a), 
X4f(T, O) imposes a soft cutoff on the energy exchange 
111 5J( d \ which appears as the last factor in Eq. {J) in 
terms of the running band cutoff D < D . In 5J^ 
(Fig. [2]b)) there is no such cutoff; however, ^/(T, fi) 
restricts the conduction electron response to the Fermi 
energy, and suppresses 5J( ex ^ compared to 8J^ by an 
overall factor of y/T^ + T 2 /Dq, as seen in Eq. The 
one-loop RG equation for the local spin exchange cou- 
pling, including RKKY corrections, Eqs. ([1]) and (JSJ) , 

reads, dJ/d\nD = -2N(0) [J + SJ^ + 5J^f '. Note 
that in this RG equation the couplings gi on sites i ^ 
are not renormalized, since this is already included in the 
full susceptibility Xif- The essential feature is that in the 
low-T limit (T <C Tk), for which Tk is to be extracted, 
the direct RKKY correction 5,j( d \ Eq. |T]), is inversely 
proportional to the renormalized Kondo scale Tr-(j/) it- 
self via %4/(0,0). Tk(v), including perturbative RKKY 
corrections, is defined as the scale where the solution of 
the RG equation diverges. This leads without further 
approximations to a highly non-linear renormalization of 
Tk given by the selfconsistency equation, 



T K (0) P \ \2g 



f(u) 
1 " /(«) 



(3) 



with g = N(0)J, f(u) =u- u 2 /2, u = yg 2 D /[±T K {y)]. 
The single-ion Kondo scale without RKKY coupling is 
Tr-(O) = Dq exp[— l/2g]. We have verified that this per- 
turbative RG treatment of the RKKY corrections is con- 
trolled in the sense that the effective perturbation param- 
eter f(u) <0.1 for all selfconsistent solutions, i.e., the ex- 
ponent in Eq. © remains small. The solution of Eq. ([3]) 
is universal in the dimensionless variables Tk{v) /Tr-(O) 
and y/y m and turns out to exist only for y smaller than 
a maximum value y m . For y/y m — ► 1 it behaves as 
TK(y)/TK (0) = r + by/l — y/y m , with numerically deter- 
mined constants r = T K (y m )/T K (0) w 0.391, b w 0.517. 
(cf. Fig. [2t))- The maximum RKKY parameter y m be- 
yond which a solution of Eq. ([3]) ceases to exist, can be 
deduced in terms of tk = Tk{0)/D o as, 



\tt k (lnrx) 2 



2 -'"t-#->°t)' 



(4) 



(Fig. [31 red line). For y>y m , the Kondo spin exchange 
coupling including RKKY corrections does not diverge 
under RG. Hence, as the essential result of this analy- 
sis, in a dense system of Kondo ions, complete Kondo 
screening ceases to exist above a critical RKKY coupling 
strength, y>y m , even if magnetic ordering does not oc- 
cur. For comparison, the well-known Doniach criterion 
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FIG. 3: (color online) Schematic T — phase diagram of a 
HF system with a magnetic QPT driven by RKKY coupling 
or a SDW instability in the plane of single-ion Kondo-scale 
Tk(0)/Do and RKKY coupling strength y, as drawn from the 
high-T analysis. Red line (ym)'- maximum RKKY coupling 
strength beyond which complete Kondo screening ceases to 
exist, Eq. Q. At the red line a step occurs in Tjj-(j/). Blue line 
(j/sdw): RKKY coupling strength at which the heavy Fermi 
liquid becomes unstable towards a SDW, assuming a spin-spin 
coupling Jsdw = J- The cases y m < j/sdw and y m > j/sdw 
distinguish whether the QPT is LQC-like (yellow) or HM-like 
(white), respectively. The arrows indicate estimates |l6j ] for 
Tk(0)/D o for CeCus-zAu^ and for Ce]SH 2 - x Cu x Ge 2 . 

[28| pertains to the breakdown of Kondo screening due 
to magnetic ordering and reads Tr-(O) r* y m N(0)J 2 . It 
is similar, albeit not equivalent, in numerical values to 
Eq. ([JJ and does not allow for a determination of the 
composition-dependent screening scale Tk(u)- By iden- 
tifying y m with the QCP, as done in Fig. [21 the x depen- 
dence of Tf[[y(x)] near this critical cndpoint explains the 
abrupt step of Tk observed experimentally in the high- 
T spectral data. The small KR spectral weight seen for 
x > x c should be interpreted merely as the logarithmic 
onset of Kondo screening which does not fully develop 
even for T — > 0. Consequently, the small, residual local 
moments must order in dimensions d > 2 at sufficiently 
low T. 

Hence, the theory predicts two possible scenarios, de- 
pending on the size of Tk (0) / Do: (I) Magnetic ordering 
at T = occurs for an RKKY parameter y = ysnw < Vmi 
i.e., without breakdown of Kondo screening. The order- 
ing may arise in this case from a T = spin-density-wave 
(SDW) instability of the underlying heavy Fermi liquid. 
This corresponds to the HM scenario, depicted in Fig. [3J 
as the white region. In this case, a step-like behavior 
of Tk as extracted from the high-T UPS spectra is still 
predicted from Eq. ([3]) at y = y m , i.e., inside the region 
where magnetic ordering occurs at T < Tv. (2) Mag- 
netic ordering does not occur for y < y m . In this case 
Eq. indicates a breakdown of Kondo screening at the 
magnetic QCP, and quantum critical fluctuations (not 
considered in our theory) will suppress the low-T spin 
screening scale to zero at this point. This is the LQC 
scenario, shown in Fig. [3] as the yellow region. 



To conclude, the theoretical analysis predicts gener- 
ally that an abrupt step of the Kondo screening scale 
extracted from high-T spectral data should occur in any 
HF compound with competing Kondo and RKKY inter- 
actions, as long as the single-ion Kondo screening scale is 
larger than the magnetic ordering temperature. Whether 
this distinct feature is located at the quantum critical 
control parameter value x c or inside the magnetically or- 
dered region constitutes a general high-T criterion to dis- 
tinguish the LQC and HM scenarios. The sharp step of 
Tk occurring in our UPS data very close to x c (see Fig. [5]) 
strongly supports that CeC\iQ- x Au x falls into this class, 
as was previously inferred indirectly from inelastic neu- 
tron scattering data [l4[ ■ 

We would like to thank F. Assaad, S. Kirchner, and 
M. Vojta for fruitful discussions. This work was sup- 
ported by DFG through Re 1469/4-3/4 (M.K., A.N., 
F.R.), SFB 608 (J.K.) and FOR 960 (H.v.L.). 



[1 
[2 

[3 

[4 
[5 

ir 

[8 
[9 

[10 

[11 
[12 

[13 
[14 
[15 
[16 
[17 
[18 
[19 

[20 

[21 
[22 
[23 
[24 
[25 



[26 
[27 
[28 



1165 (1976). 
7183 (1993). 

Smith, Nature 413, 804 



H. v. Lohneysen, A. Rosch, M. Vojta, and P. Wolfle, Rev. 
Mod. Phys. 79, 1015 (2007). 

P. Gegenwart, Q. Si, and F. Steglich, Nature Phys. 4, 
186 (2008). 

J. A. Hertz, Phys. Rev. B 14, 

A. J. Millis, Phys. Rev. B 48, 
Q. Si, S. Rabello, and J. L. 
(2001). 

P. Coleman, C. Pepin, Q. Si, and R. Ramazashvili, J. 

Phys.: Condens. Matter 13, R723 (2001). 

H. v. Lohneysen et al., Phys. Rev. Lett. 72,3262 (1994). 

H. v. Lohneysen, J. Phys. Cond. Mat. 8, 9689 (1996). 

H. v. Lohneysen et al, J. Magn. Mag. Mat. 177-181, 12 

(1998). 

O. Stockert et al, Phys. Rev. Lett. 80,5627 (1998). 
H. v. Lohneysen et al, Eur. Phys. J. B 5, 447 (1998). 

B. Stroka et al, Z. Phys. B 90, 155 (1993). 

K. Grube et al, Phys. Rev. B 60, 11947 (1999). 

A. Schroder et al, Nature 407, 351 (2000). 

F. Reinert et al, Phys. Rev. Lett. 87, 106401 (2001). 

D. Ehm et al, Phys. Rev. B 76, 045117 (2007). 

M. Klein et al, Phys. Rev. Lett. 101, 046406 (2008). 

J. Kroha et al, Physica E 18, 69 (2003). 

H. v. Lohneysen et al, J. Alloys Compd. 408—412, 9 

(2006). 

H. G. Schlager, A. Schroder, M. Welsch, and 
H. v. Lohneysen, J. Low Temp. Phys. 90, 181 (1993). 
F. Patthey et al, Phys. Rev. B 42, 8864 (1990). 
M. Gamier et al, Phys. Rev. Lett. 78, 4127 (1997). 
J. Allen et al, et al., J. Appl. Phys. 87, 6088 (2000). 
D. Ehm et al, Acta Phys. Pol. B 34, 951 (2003). 
T. Costi, J. Kroha, and P. Wolfle, Phys. Rev. B 53, 1850 
(1996); N. E. Bickers, D. L. Cox, and J. W. Wilkins, 
Phys. Rev. B 36, 2036 (1987). 

B. A. Jones and C. M. Varma, Phys. Rev. Lett. 58, 843 
(1987). 

H. Maebashi, K. Miyake, and C. M. Varma, Phys. Rev. 

Lett. 95, 207207 (2005). 

S. Doniach, Physica B & C 91, 231 (1977). 



